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Abstract 

Applying the recently developed formalism of quantum-state transforma- 
tion at absorbing dielectric four-port devices [L. Knoll, S. Scheel, E. Schmidt, 
D.-G. Welsch, and A.V. Chizhov, Phys. Rev. A 59, 4716 (1999)], we calcu- 
£SJ . late the quantum state of the outgoing modes of a two-mode squeezed vacuum 

transmitted through optical fibers of given extinction coefficients. Using the 
Peres-Horodecki separability criterion for continuous variable systems [R. Si- 
mon, Phys. Rev. Lett. 84, 2726 (2000)], we compute the maximal length of 
transmission of a two-mode squeezed vacuum through an absorbing system 
& for which the transmitted state is still inseparable. Further, we calculate the 

maximal gain for which inseparability can be observed in an amplifying setup. 
Finally, we estimate an upper bound of the entanglement preserved after trans- 
mission through an absorbing system. The results show that the characteris- 
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tic length of entanglement degradation drastically decreases with increasing 
strength of squeezing. 



X 

1 Introduction 



Recently there has been increasing interest in the use of entangled continuous variable 
systems in quantum communication |l] . The most prominent example is the two- mode 
squeezed vacuum which in the Fock basis reads as 

CO 

|^ = e c(44-aia 2 )| 00 ) = ^/YZ^J^q" 1 \nn) (1) 

n=0 

[q = tanh^, ( real] and whose Wigner function is a Gaussian. When transmitted 
through a noisy channel, such as optical fibers of given extinction coefficients, the 
resulting state will in general be mixed. For further use of the state in some quantum 
communication experiment the question of the entanglement that is available after 
the transmission arises. Real entanglement measures, however, are difficult to com- 
pute and are mostly known only numerically. A typical example is the entanglement 
measure based on the relative entropy where the distance of the density matrix of the 
state under consideration to the set of all separable density matrices must be com- 
puted [[|. In practice, this becomes impossible for states like a two-mode squeezed 
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vacuum with reasonable strength of squeezing, because for comparable numerical ac- 
curacy the Hilbert space can only be truncated at higher dimension with increasing 
squeezing strength. 

Here, we proceed in a different direction. Starting from the quantum-optical input- 
output relations of light at absorbing dielectric four-port devices Q, we first present 
the basic formulas for determining the output quantum state from the input quan- 
tum state and the characteristic transmission and absorption matrices of the devices 
(Sec. 0). We then apply the input-output formalism to the calculation of the output 
Wigner function observed when two modes that are initially prepared in a squeezed 
vacuum are transmitted through absorbing fibers. Using the Peres-Horodecki separa- 
bility criterion for Gaussian states we calculate the maximal length of transmission 
for which the output state is still inseparable in principle (Sec. §). Calculating the 
density matrix of the transmitted light in the Fock basis, on applying the formalism 
developed in p, we finally extract some pure state from the output density matrix, 
use the convexity property of the entanglement measure 0, and derive an estimate 
of the amount of entanglement available after transmission (Sec. 



2 Quantum- state transformation 

Let us first consider the quantum-optical input-output relations and the correspond- 
ing quantum-state transformation formulas for light at dispersing and absorbing di- 
electric four-port devices. We restrict ourselves to a quasi one-dimensional scheme, as 



Figure 1: Scheme of the dielectric four-port device, with di = ai(uj) 
[bi = t>i(u)] and §i = gi(uj) being the destruction operators 
of the input [output] photons and the device excitations 
respectively. 

depicted in Fig. [I], in which the dielectric device is surrounded by vacuum. Applying 
the formalism developed in [0, we quantize the electromagnetic field in the presence 
of the device by means of a Green function representation of the field and introduction 
of bosonic fields playing the role of the collective excitations of the field, the dielec- 
tric matter, and the reservoir. It turns out that outside the device the usual mode 
expansion applies, with the CLiioS) and bi(u) in Fig. [l| being respectively the photonic 
operators of the incoming and outgoing plane waves at frequency u. The cji{uj) in the 
figure are the bosonic operators of the device excitations and play the role of noise 
forces associated which absorption. It then follows that the action of the dielectric 
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device on the incoming radiation can be described by quantum optical input-output 
relations which, in fact, are nothing but a suitable rewriting of the corresponding 
one-dimensional Green function. Let us introduce the compact two-vector notation 
a(uj), h(u), and g(u;) for the field and device operators respectively. The input-output 
relations can then be written in the compact form || 

b(w) = T(w)a(w) + A(u)g(u), (2) 

where the characteristic transmission and absorption matrices T(u) and A(u;), re- 
spectively, satisfy the energy- conservation relation 

T(co>)T + (u;) + A(uj)A + (uj) = I. (3) 

Equations (0) and (|3]) are valid for any frequency. They allows us to construct the 
output operators from the input operators over the whole frequency range. 

The second step in the quantum-state transformation corresponds to an open- 
systems approach. Suppose the incoming field is prepared in some state of the Hilbert 
space H&eid and the device (including the reservoir) is initially prepared in some state 
of the Hilbert space Hdevice- In the full Hilbert space, which is the tensor product 
Afield <8> Tidevice, a unitary operator transformation can then be constructed, whereas 
in the space Hfi e id it could not due to the dissipation processes. Let us define the 
four-vector operators 



Ct(UJ) 



SL{UJ 



where h(u) is some auxiliary (two-vector) bosonic device operator. Then, the input- 
output relation (§) can be extended to a four-dimensional transformation 

= A(w)d(w) , (5) 

with A{uj) e SU(4) [|]. Explicitly, 

A(u)-( TM AM (6) 

{U) \ -S{u)C-\u)T(u) C(u)S-\u)A(u) J W 

with the commuting positive Hermitian matrices 



C(u) = y/T(u)T+(u), S(u) = y/A(ujA~Hu) . (7) 
Hence, there is a unitary operator transformation 

j9(w) = ffla(u))U (8) 

where 

?7 = expj-z Jdu [a f (u;)] T *(u;)d(u/)J (9) 
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and 

A(w) = e^^ . (10) 

Note that U acts in the product space Hadd <8> Tidevice- Given a density operator o in of 
the input quantum state as a functional of 6t(u), the density operator of the output 
quantum state is obtained by a unitary transformation with the operator U from 
Eq. (H) and projecting back onto the Hilbert space 7ifi e id- Hence, 

= Tr (D) {f/fe^} = TrW {£ in [A + ( W )d( W ), A»aV)] } , (11) 

where Tr^ - 1 means tracing with respect to the device variables. Note, that the dif- 
ference to usually considered open-systems theories is provided by the fact that we 
actually know how the dissipative environment (e.g., a dispersing and absorbing fiber) 
acts on our quantum states. 

Let us briefly comment on amplifying devices. In contrast to absorbing devices, we 
have now to insert the noise creation operators gj(uj) into the input-output relation 
dD [i.e., g(u) — >gt(a>)]. The relation (|3|) then changes to 

T(u)T+(u) - A(w) A + (w) = I , (12) 

where A(u) plays the role of the gain matrix. Further, the 4x4-matrix A(u) becomes 
an element of the noncompact group SU(2,2). 

3 Application of the Peres— Horodecki separability 
criterion 

Let us consider a two-mode squeezed vacuum which is transmitted through a noisy 
communication channel such as two absorbing (amplifying) dielectric four-port de- 
vices as sketched in Fig. Q the characteristic transmission and absorption (gain) 
matrices of the devices being T™(w) and A™(oj) respectively (i — 1,2). In Eq. ([!]), 
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Figure 2: A two-mode input field prepared in the quantum state 
\ip) is transmitted through two absorbing (amplifying) 
dielectric four-port devices, Si, S3 {a' 2) a'^ being the pho- 
tonic operators of the relevant input (output) modes. 

the two-mode squeezed vacuum is given in the Fock basis. Equivalently, it can be 
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expressed in terms of the Gaussian Wigner function 

W{£) = (47rVdetv)~ exp-j-f^V^^} 



(13) 



where £ is a four-vector whose elements are the quadrature-component variables 
Qi,Pi, Q3,P3, and V is the 4x4 variance matrix of the Wigner function, 



V 



X 



z 

Y 



/ c/2 -s/2 \ 

c/2 s/2 

-s/2 c/2 

\ s/2 c/2 / 



(14) 



(c = cosh2C, s = sinh2£). Transmitting the two- mode squeezed vacuum through the 
four-port devices at some temperatures the Wigner function of the transformed 
state is again a Gaussian. Using the input-output relations for absorbing (amplifying) 
devices as given in Sec. |2], we can easily transform the input variance matrix (jnj) to 
obtain the output variance matrix. The result is 



\^a! 2 ) s 
(a' 2 a' 4 ) 



WlS/At; 



T$\ 2 {a%) 8 +\T£\\a\a±) 



(2)|2/ A t; 



(2),2/W2)t^(2)\ 



(15) 



(the subscript s refers to symmetric operator ordering as required for the Wigner 
function) which can also be rewritten as correlations of the quadratures. Introducing 
the abbreviating notation =Tj and T 22 =Ri, we arrive at the following elements 
of the output variance matrix: 

X n = X 22 = \c\Tr\ 2 + \\R^ + a (n thl + \) (l - |T\| 2 - {R^ 2 ) 

\c\T 2 \ 2 + \\R 2 \ 2 + a (n th2 + ±) (l - \T 2 \ 2 



Y u = Y 



22 



\R2\ 2 ) 



-± S Re (TxT 2 ), 
-±slm (TxT 2 ) 



(16) 

(17) 
(18) 
(19) 



■Zll — — -^22 
Z\2 = -^21 r 

(X 1 2 = X 2 i = Y 1 2 = 5 / 2i = 0), where cr = +l(— 1) for absorbing (amplifying) devices, and 
^th« = [exp huj/{kB'&i) — l]^ 1 is the mean number of the (thermal) excitations of the 
ith device. 

Application of the Peres-Horodecki separability criterion j| 

, 2 



detXdetY+ (± - | det Z|) - Tr (XJZJYJZ T j) > ±(detX + detY) 



with 



1 
-1 



(20) 
(21) 
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to the output variance matrix yields (for equal devices) the inequality 

^> ^^mfc^ . (22, 



Hence, for chosen squeezing parameter ( and chosen transmission and reflection coeffi- 
cients T and R, respectively, there exists a maximal temperature and correspondingly 
a maximal mean excitation number of the thermal state in which each of the two de- 
vices is prepared such that the quantum state of the transmitted squeezed vacuum 
is still not separable. 

In particular, for absorbing devices (cr = +l) the inequality (|22|) reads 





T 


2 (1" 


- e-^l) 


2 


(1 




R 


2 _ 


T 


2 ) 



"th > ' ' V |D| , , ' • (23) 



With regard to optical fibers with perfect input coupling, we may let R = and relate, 
according to the Lambert-Beer law, the transmission coefficient to the propagation 
length I as \T\ = e~ l/lA , with l A being the characteristic absorption length of the fibers. 
From the inequality (^) it then follows that the upper bound of the propagation 
length for which the transmitted squeezed vacuum is still not separable is 



I = hi a In 



1 + — (1 -e- 2 KI) 
2n th y > 



(24) 



It is interesting to note that this result agrees with that calculated in ||, if the 
'renormalized time' in |S| is replaced with 1 — |T| 2 . Note that in our formalism the 
bound is simply obtained by using the quantum-optical input-output relations (0) for 
light at dispersing and absorbing dielectric four-port devices. Involved calculations of 
the dynamics of the Wigner function are not needed. 

Moreover, the general result (^) also applies to amplifying devices (cr = — 1). In 
particular, in the zero-temperature limit (n t h = 0) the boundary between inseparability 
and separability is reached when 

, l9 2(l-|i?| 2 ) , . 

T = rz -2i i 25 

1 + e ^ 

is valid. For zero reflection (i? = 0), Eq. fl25|) reveals that the upper limit of the 'excess' 
gain g = |T| 2 — 1 > for which inseparability changes to separability is simply given 
by the squeezing parameter q, 

g= \q\ =taoh|C|. (26) 

An obvious consequence of Eq. ( pip is that entanglement cannot be produced from 
the vacuum by amplification. For the vacuum the squeezing parameter has to be set 
equal to zero and thus from Eq. (|26| ) it follows that any nonvanishing gain g must 
necessarily lead to a separable state. Another interesting fact is that there exists an 
absolute upper bound of the gain for which inseparability can be retained. Since the 
absolute value of the squeezing parameter (and thus the 'excess' gain g) is bounded 
by 1, one is left with the conclusion that an amplifier which doubles the intensity of a 
signal (|T| 2 = 2) destroys all but the maximal entanglement in a two-mode squeezed 
vacuum state which spoils the use of fiber amplifiers in quantum communication. 
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4 Entanglement estimates 

The separability criterion exploited in Sec. [I] can tell us only if the transmitted quan- 
tum state is separable or not. It can not, however, provide us information about the 
amount of entanglement which is actually contained in the state. In order to obtain 
analytical estimates of the available amount of entanglement, we note that any en- 
tanglement measure (such as the distance of the state under consideration to the set 
of separable states measured by the relative entropy) has the convexity property || 

£7[A£i+(l-A)&] < XEfa) + (1 - X)E{q 2 ). (27) 

This property can advantageously be used to find bounds on the entanglement. If 
we are able to divide the quantum state into a sum of separable states (having no 
entanglement) and a single pure state, then an upper bound on the entanglement is 
given by the reduced von Neumann entropy of the extracted pure state |J. 

Let us consider two modes that are initially prepared in a truncated version of 
the quantum state (|l|) 

|0( ff )> =i — = (I00> H- g|ll» , (28) 

l + q 2 



which approximates a two-mode squeezed vacuum for small values of the squeezing 
parameter (i.e., \q\ <1). It is not difficult to prove that the entanglement of the state 
is 

E = Ul + q 2 ) -zr—5 In q 2 ■ (29) 
1 + q l 

Applying the transformation formula (]TT|) , the quantum state in which the two modes 
are prepared after propagating through two fibers of transmission coefficients T\ and 
T 2 is derived to be 

£1 = [(1 " \Ti\ 2 ) (1 - |T 2 | 2 ) |00)(00| + \T,\ 2 (1 - |T 2 | 2 ) |10)(10| 



where 



1 + 1 


ft 




1 + 


lei 


2 



+ |T 2 | 2 (1 - \T,\ 2 ) |01)(01|] + W))W% (30) 



q' = T{T 2 q. (31) 

Here and the following we assume that the fibers are prepared in the ground state 
(low-temperature limit). Since the first term on the right-hand side in Eq. ( p0[ ) is a 
sum of separable states and the second term is a pure state whose entanglement is 
given by Eq. ( [29]) with q' in place of q, the entanglement of the state in Eq. ( |30"|) can 
be estimated, on recalling the convexity property (p7|), according to 

E < [(l + mm 2 ) In (1 + \qT,T 2 \ 2 ) - \qT x T^\n \qT,T 2 \ 2 ] . (32) 
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We see that with increasing propagation lengths, i.e., with decreasing transmission 
coefficients, the entanglement of the transmitted light decreases more rapidly than 
the intensity. 

Now let us return to the exact two-mode squeezed vacuum. Applying Eq. (|TTD 
and transforming the density operator g m = with from Eq. (|]), we derive 

after a lengthy, but straightforward calculation 



oo m 1 



m=0 n=0 



mini 



5^Cm,n,fc (c m _ n \m-n + k)(k\ + E.c.) 



k=0 



^Drn,n,l (d m - n \m ~ U + /) (l\ + H.C.) 



1=0 



where 



C m ,n,k 



q n m\n\ (1 — \T\ 



2k 



(n — k\)yjk\{m — n — k\) 

(m-n)/2rpm-n fi _ If \ 

y ± i \ L 2 mn ) 



d„ 



q n m\n\ (1 — \T 2 



To 



12/ 



(n-l\)^/l\(m-n-l\) 
= q (™-">V*T™- n (1 - \5 mn ) 
Performing the sum over n then yields the final result 



oo oo 



m=0 k,l=0 



where 

^■k,l,m 



(33) 

(34) 
(35) 
(36) 
(37) 

(38) 



[g 2 (l-|Ti| 2 )(l-|T 2 | 2 )]°a!(a+m)! / \T X 



^kll\(k+m)\(l +m)\(a — k)\(a — l)\ 







2 


1- 




T\ 


2 





T 2 


2 


1- 




T 2 





X 3^2 



a + I, a + m + 1,1 



(39) 



[a = max(fc, I)]. Note that one lower index of the hypergeometric function 3F2 is always 
equal to unity, so that we effectively deal with a Gaussian hypergeometric function 

2^1- 

We now try to decompose the density operator fl3"8|) into a pure state \^) and 
some residual state g' whose entanglement is desired to be sufficiently small, 

eS = A^ + (i-A)|*><*|, 



A suitable pure state |^) may be chosen such that 



y/l - A|*> 



000 



^ K 00n c n d n \nn) . 



(40) 



(41) 



n=0 
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It has the properties that (i) only matrix elements of the same type as in the initial 
squeezed vacuum occur and (ii) the coefficients of the matrix elements 1 00) <-> \nn) 
are met exactly, i.e., 

(1 - A)(00|*)(*|nn) = (00\g<£>\nn). (42) 

This choice is of course not unique. Moreover, the residual state might still contain 
some entanglement. For small enough squeezing parameter q, however, the residual 
entanglement is expected to be small compared to the entanglement contained in the 
state i.e., XE(g') <C (1 — X)E(\^f)). In principle, one can proceed and extract 

more and more pure states from the output state and apply the generalized inequality 



E 



Disregarding a possible (small) entanglement of the residual state g', the en- 
tanglement of the pure state |\?) gives us some estimate of (the upper bound of) 
entanglement of the output state g ou [, 



E(gS t ) * (1 - X) E(\*)) 

1 —x 



1 — x 

In 



•a;) 2 - 



[l-x) 2 -y 



(l-x){[y + (l-x) 2 ] \n(l-x)-y\ny} 

[y-{l-x) 2 f ' 1 > 



where 



x = g 2 (l-|T 1 | 2 )(l-|T 2 | 2 ), (45) 
y=\qT 1 T 2 \ 2 . (46) 
Note that for T\ = T 2 = 1 Eq. ([0]) gives the correct entanglement of the input state, 



E(^t) = m)) = - In (1 - q 2 ) - j-^r In q\ (47) 

Ti=t 2 =i (1 — q z ) 

In Fig. |3] we have illustrated the dependence of the estimated entanglement of 
the transmitted two- mode squeezed vacuum, Eq. (0), on both the fiber length and 
the strength of the initial squeezing. It is seen that with increasing strength of the 
initial squeezing [which is, according to Eq. (|47|) , a measure of the strength of the 
initial entanglement] the entanglement of the transmitted light drastically decreases 
with the transmission length. The transmission length at which the entanglement 
degradation has reached half of the initial entanglement is shown in Fig. [| as a 
function of the mean number of initial photons (n) , which is related to the squeezing 
parameter \q\ 2 according to 



2 



> (n) 

q= WTi- (48) 

Figure (| reveals that even for relatively small mean number of photons this charac- 
teristic length of entanglement degradation is much more shorter than the absorption 
length. Obviously, entanglement cannot be maintained when going to more macro- 
scopic nonclassical states. Since strong squeezing (i.e., large photon number) is typ- 
ically required in quantum teleportation |], one has to find a compromise between 
highest possible entanglement and lowest entanglement degradation. 
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Figure 3: Estimate of the entanglement of the transmitted squeezed 
vacuum, Eq. (33), as a function of the initial squeezing 
parameter \q\ 2 and the transmission length I [E(\ip)), ini- 
tial entanglement]. 




10 15 20 25 30 35 



40 



(n) 



Figure 4: Entanglement degradation length Ze after which the en- 
tanglement of the transmitted two-mode squeezed vac- 
uum has decreased to half the initial value as a function 
of the initial value of the mean number of photons (n) 
(I a, absorption length). 



5 Summary and Outlook 

We have studied the entanglement degradation of a two-mode squeezed vacuum state 
transmitted through noisy communication channels, a typical example being the prop- 
agation along dispersing and absorbing or amplifying optical fibers. As expected, both 
absorption and amplification lead to entanglement degradation, because of the addi- 
tional noise introduced by them. Using the quantum-optical input-output relations 
of radiation at dielectric four-port devices, we have derived the maximal transmis- 
sion length after which an initially entangled state becomes separable. Analogously, 
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we have found that there is a maximal gain factor for which inseparability (of any 
two-mode squeezed vacuum state) can be retained at all. 

Knowledge of the complete quantum state at the output of the device enables 
us, in principle, to compute the amount of entanglement available after transmission 
such as the distance of the state to the set of all separable states. However, for 
higher- dimensional Hilbert spaces as in our case this is impossible to do in general. 
Therefore we have restricted ourselves to the calculation of upper bounds or estimates 
of upper bounds. For weak squeezing we may truncate the Hilbert space effectively 
at low photon numbers providing us with a way to establish an upper bound on the 
entanglement by exploiting the convexity property of entanglement measures. The 
procedure is based on the extraction of a single pure state from the output state 
leaving behind only separable states. For the general case of infinite dimensional 
Hilbert space we have derived an estimate of the upper bound. 

Still, we are left with (estimates) of bounds on the entanglement. Future works 
must surely contain algorithms for computation of the entanglement, at least for 
Gaussian states. A possible step in that direction would include the calculation of 
the distance of a given Gaussian state to the set of all separable Gaussian states whose 
surface can again be parametrized by the Peres-Horodecki separability criterion. 
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